TYPE lib SUPERNOVAE WITH COMPACT AND 
EXTENDED PROGENITORS 
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ABSTRACT 

The classic example of a Type lib supernova is SN 1993J, which had a cool 
extended progenitor surrounded by a dense wind. There is evidence for another 
category of Type lib supernova which has a more compact progenitor with a 
lower density, probably fast, wind. Distinguishing features of the compact cate- 
gory are: weak optical emission from the shock heated envelope at early times; 
nonexistent or very weak H emission in the late nebular phase; rapidly evolving 
radio emission; rapid expansion of the radio shell; and expected nonthermal as 
opposed to thermal X-ray emission. Type lib supernovae that have one or more 
of these features include SNe 1996cb, 2001ig, 2003bg, 2008ax, and 2008bo. All 
of these with sufficient radio data (the last four) show evidence for presupernova 
wind variability. We estimate a progenitor envelope radius ~lx 10 11 cm for SN 
2008ax, a value consistent with a compact Wolf-Rayet progenitor. Supernovae in 
the SN 1993J extended category include SN 2001gd and probably the Cas A su- 
pernova. We suggest that the compact Type lib events be designated Type cllb 
and the extended ones Type ellb. The H envelope mass dividing these categories 
is ~ 0.1 M Q . 

Subject headings: Circumstellar matter — shock waves — supernovae: general 
- supernovae: individual (SN 1993J) — supernovae: individual (SN 2008ax) 

1. INTRODUCTION 

A Type lib supernova has a spectrum with high velocity H lines near maximum, but in 
the late nebular phases more closely resembles a Type lb supernova, with weak or absent H 
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lines. The first supernova placed in this category was SN 1987K (jFilippenkolll988l ) . However, 
it was the nearby SN 1993J in M81 that has come to best define the properties of a SN lib . 
In this case, the progenitor star was identified as a K supergiant (lAldering et al.l Il994j ). 
consistent wit h the interpretation o f the early light curve which required a stellar radius of 
~4x 10 13 cm ( IWoosley et al.lll994l ). Such a star is expected to have a slow, dense wind, and 
there is radio and X-ray evidence for a wind wit h mass loss rate M^4x 10~ 5 yr -1 for 
an as sumed wind velocity of v w = 10 km s _1 (IFransson et al.lll996l ; iFransson fc Bjornsson 



9981). Radio VLBI observations provi de strong support for the wind interaction picture 



(IMarcaide et al.lll997l ; iBartel et al.ll2002l ). Th e interaction with the dense wind also produces 



some Ha emission in the late nebular phase (IMatheson et al.ll2000l ) 



A consistent scenario for SN 1993J is thus the explosion of an extended star that interacts 
with the dense wind from the progenitor star. However, there is also evidence for Type 
lib supernovae arising from compact stars, and that is th e topic addressed here. Some 
of that evidence has al ready been present ed for SN 200 lig (IRvder et al.l 120041 ) . SN 2003bg 
JSoderberg et all l2006h . and SN 2008ax rtPastorello et al.l l2008h . In § 2. we present the 
differences between SNe lib that have extended or compact progenitors. The results are 
discussed in § 3. 



2. TYPE lib SUPERNOVA CATEGORIES 



2.1. The Progenitor Star and Early Light Curve Evolution 



The most direct way to observe the progenitor star is in presupernova observations. 
This has bee n possible for SN 1993 J, as discussed above, and has recently been possible for 
SN 2008ax. ICrockett et aD J2OO8I ) find that the progenitor of the Type lib SN 2008ax is 
either a single, He-rich, Wolf-Rayet star or a stripped star in an interacting binary in a low 
mass cluster. 

The light curve evolution in the first days after the explosion, which is the time of cooling 
envelope emission after shock breakout, provides a good indicator of the initial stellar extent. 
The cooling envelope emission is more luminous for a more extended star because of the larger 
emitting area and adiabatic expansion losses are less important for an extended star. The 
early peak luminosity and first 10 days of ev olution provided evid ence for a large radius, 
~4x 10 13 cm for the progenitor of SN 1993J (IWoosley et al.lll994l ). 



In the case of SN 2008a x, a comparable early luminous phase was absent in the light 
curve (jPastorello et al.l 120081 . Fig. 11). Very early o bservations with Sw ift did show evi- 
dence for an initial peak in the uvwl, u and h bands ( jRoming et al.l 120091 ). The luminosity 
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level of the emission is comparable to, although slightly fa inter than, similar emission ob- 
served in the Type Ib/c SN 199 9 ex (jStritzinger et al.l 120021 ) . We used the cooling envelope 
theory of I Chevalier fc Franssonl (120081 ) to model the uv w2, uvwl, u, b emission on day 1.7 



jRoming et alJl2009h and the q, r l i,z e mission on day 1.9 (IPastorello et ajj|2008j ). The mag- 



nitude results of iRoming et al.l (120091 ) were converted to fluxes using the flu x conversion 



factor s, pivot wavelengths, and bandpass widths for the UVOT as reported in iPoole et al. 



(120081 ). For a distance of 9.6 Mpc, a blackbody fit to the fluxes gives an effective temperature 
T ~ 6000 K and a photospheric radius r p h = 4.3 x 10 14 cm on day 1.8. Taking an ejec t a mas s 
of M e j = 3.8 M q and an energy E = 1 . 5 x 10 51 ergs from the model of iTsvetkov et al.l (120091 ). 
equation (2) of Chevalier fc Fransson ( 20081) yields r p h. = 3.8 x 1 14 cm , in good agreement 



with the observed value. Equation (5) of I Chevalier &: Franssonl (120081 ) and the observed T 
yield a progenitor radius of R = 1 x 10 11 cm. However, the t emperature is sufficiently low 
that recombination is becoming important, so that the model of lChevalier &: Franssonl (120081 ) 
is inaccurate. The progenitor radius estimate should still be valid to within a factor 2 — 3, 
showing that the progenitor was compact, with a radius consistent with a Wolf-Rayet star. 
The radius determined here is that of the stellar envelop e; the ph otospheric radius may be 
a factor 2 — 3 larger due to the stellar wind (see Fig. 1 of lLill2007l). Althou g h the explosion 



param eters for SN 2008ax are similar to SN 1993J (IPastorello et al.l 120081 : ITsvetkov et al. 



20091 ). the progenitor is very different, as indicated by the progenitor observations and the 
early light curve. 

For a compact star, the cooling envelope emission lasts for only a few days, so that 
observations must be carried out very soon after the explosion for this phase to be detected. 
Observations that rule out an extended progenitor may be the most that can be obtained. For 
the Type lib SN 1996cb, there are light curv e observations s howing that an early luminous 
phase like that of SN 1993 J was not present (IQiu et al.lll9991 ). 



2.2. Circumstellar Interaction 



Circumstellar interaction depends on the progenitor radius in that more compact stars 
have faster, and thus more rarefied, winds and shock acceleration at the time of breakout 
extends to higher velocities in compact stars. Radio emission is an excellent indicator of cir- 
cumstellar interaction and Fig. 1 shows peak radio luminosity vs. time of peak for SN lib and 
SN Ib/c. SSA (synchrotron self- absorption) is expected to be the dominant absorption mech- 
anism for the high velocities present in SN Ib/c a nd can be used t o estimate the radio shell 
velocity at the time of the peak radio lum inosity (IChevalierlll998f) . The SN l ib designation 
has been given to SN 1993J, SN 2001gd jMatheson et al.lboOlh . SN 2001ig Jphillips et al. 
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20011), SN 2003bg (Irlamuv et al.l 120091 ) . and SN 2008ax flPastorello et al.l l2008h . Fig. 1 



shows that there is a class of SN lib that have radio properties similar to those of SN 
Ib/c, i.e. they have rapid radio evolution for a given luminosity, indicating velocities of 
30, 000 — 50, 000 km s _1 . These velocities are expected for the explosion of a Wolf-Rayet star 
because of strong shock acceleration in the outer parts of the sta r and the low wind density 
resulting from a high wind velocity (j Chevalier fc Franssonl 120061 ). They are designated here 
as SN cllb because of the expectation that they have relatively compact progenitors. There 
is direct light curve and progenitor evidence for this in the case of SN 2008ax, as discussed 
above. SN 1993J is placed in the SN ellb category because of its extended progenitor, as is 
SN 2001gd, be cause of its radio prop erties. In addition to the radio light curve information 
on SN 2001gd (jStockdale et al.ll2007l ). there are also VLBI observations on days 286 and 582 
from explosi on, which combine to giv e a mean velocity ~ 9100 km s -1 on day 434 with 20% 
uncertainty (IPerez- Torres et al.ll2005l ). This velocity is consistent with that found from the 
radio light curves (Fig. 1) and is somewhat less than th at observed for the i nteraction shell 
in SN 1993J at a comparable time, ~ 14, 000 km s _1 (IBartel et al.l 120021 : iMarcaide et al. 



1993). 



VLBI observations of SN 2008ax dMarti-Vidal et al.ll2009l l indicate a velocity ~ 52, 000 km s" 1 
(in a plausible model) at an age of 33 days, supporting the SSA model. The h ighest H ve- 
locity in optical lines was only about half this, as noted by iMartf-Vidal et al.l (120091 ) . This 
is not a problem for the Type cllb case because the high velocity shocked H is nonradiative, 
as discussed below. The unshocked gas is of such low density t hat it is not seen. T here do 
not seem to be sufficient early radio observations of SN 2008ax (jRoming et al.l 120091 ) to test 
for SSA from the light curves , but in the case of SN 2003bg, there is support for the SSA 
model (jSoderberg et al.l 120061 ). 



The SSA model of I Chevalier fc Franssonl (120061 ) can be used to obtain estimates of 
the mass loss properties of the compact supernovae. We characterize the mass loss density 
density by A*, where M/v w = A* x 10" 5 M© yr _1 /1000 km s" 1 and find A*e B -ia 8 / 19 = 3 
for SN 2001ig, 13 for SN 2003bg, and 0.65 for SN 2008ax; here, e B _i is the fraction of the 
postshock energy density in magnetic field in units of 0.1 and a is the ratio of relativistic 
electron energy de nsity to magnetic ener g y den sity. The results for SN 2001ig and SN 2003bg 
ar e the same as i n lCheyalier &: Franssonl (120061 ) and that for SN 2008ax uses the observations 
of Roming et al. ( 20091 ). Assuming that the efficiency of production of the radio emission is 
comparable for the different objects, the variation in the radio luminosities of the supernovae 
translates to a variation in the progenitor mass loss densities, by a factor of 20. 



Other Type lib supernovae have been detected in the radio, but are not included in Fig. 
1 because the available data do not allow the maximum luminosity to be accurately deter- 
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mined. However, there is some information. SN 2008bo showed optica l properties similar 
to SN 2008ax, including deep Ha absorption ( INavasardyan et al.l 120081 ) . There were early 
radio detections of the supernova and a later (age ~ 20 days) rebrightening ( IStockdale et al. 
20081 ). At a distance of 21 Mpc, the radio luminosity was ~ 1.5 times lower than that of 



SN 2008ax and the emission became optically thin somewhat earlier, so this object extends 
the apparent wind density to slightly lower values. There are also 2 radio observations of 
the Type lib SN 19 96cb that indicate a r = 1 frequency of ~ 11 GHz on Dec. 21, 1996 
( Ivan Dyk et al.lll996l ). The observations indicate that it also lies on the L p — t p plot close to 
SN 2008ax. 



X-ray emission is another signature of circumstellar interaction. IChevalier fc Fransson 



( 120061 ) find that the X-ray emission from SNe Ib/c is probably due to a nonthermal mech- 
anism: inverse Compton emission near optical maximum light and synchrotron emission at 
later times. The general expectation is that Type cllb events should have nonthermal X-ray 
emission, while Type ellb events similar to SN 1993J should have dominant thermal X-ray 
emission, as in SN 1993J ( IFransson et al.lll996l ). There are not yet suf ficiently well observed 
Type cllb events to test this from the spectra themselves. However, iRoming et al.l (120091 ) 



found a decrease in the X-ray luminosity of SN 2008ax from 6.0 ±1.9 x 10 erg s in the 



first month to 1.4 ±0.9 x 10 erg s in the second month. This is the magnitude of decline 
that would be expected in going from an inverse Com pton component near optical maximum 



to a synchrotron component a month later (Fig. 1 of IChevalier fc Franssonll2006l ). 



Pt 

2! 



Finally, another signature of a dense circumstellar medium is optical emission from 
gas h eated and ionized by radiation from reverse shock heated gas ( IChevalier fc Fransson 
1994T ) . The optical emission can be either from freely expanding ejecta or from a dense 
shell formed as a result of radiative cooling. SN 1993J did show evidence for persistent 
Ha emission to an age of 2500 days ( iMatheson et al.l 120001 ) : the box- like line profile and 
other characteristics were consistent with circumstellar interaction. The early Ha emission 
from SN 1993J showed a decline consistent with a radioactive power source; the late steady 
phase started at an age of 300 — 350 days dHouck fc Fransso n 1996). The estimated wind 
density in SN 1993 J corresponds to A* = 400 ( IFransson et al.l Il996l ). which is a factor 
> 30 larger than the wind densities deduced here for the SNe cllb. Cooling is primarily by 
bremsstrahlung and the reverse shock is nonradiative over the times of observations for SN 
cllb. The low density and large shock radius imply that H a from circumste l lar int eraction 
should be undetectable in the late nebular phase. However, IStritzinger et al.l ( 120091 ) suggest 
that circumstellar interaction powers Ha emission in the nebular phase from the Type lb SN 
2007Y, even though upper limits on the radio emission yield A* < 0.1. The explanation is 
probably that the observations are on days 230 and 270, when radioactivity is still a plausible 
power source. Observations of Ha at times > 1 yr are needed to distinguish between the 
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steady emission that can occur in SNe ellb and the continued drop expected in SNe cllb. 



2.3. Structure in Wind 



Both SN 2001ig (IRvder et al.l 12004 ) and SN 2003bg (ISoderberg et al.l l2006[ ) showed 
structu re in their radio lig ht curves that was interpreted as density structure in the wind 
region. iRyder et al.l ( 12004 ) suggested that the structure was due to a binary co mpanion, as 



has be en observed in the pinwheel dust structure around some Wolf-Rayet stars. IRyder et al. 



( 120061 ) in fact found evidence for a possible companion star at the site of the supernova. 
However, the star is of type late-B to late-F; such a star would have too wea k a wind to have 



much effect on the strong wind deduced for the progenitor of SN 200 lig. ISoderberg et al. 



( 120061 ) argued for stellar variability rather than a binary companion, based on the similarity 
between the radio evolution of SN 200 lig and SN 2003bg. 

The additional cases of Type cllb supernovae also show evidence for wind density varia- 
tions. The radio light curv es of SN 2008ax show a peak in the optically thin phase at t ~ 60 
days ( Roming et al.l 12009) . In the case of SN 2008bo, there is a possible re-brightening at 
t ~ 20 days (jStockdale et al.ll2008l ). There is thus evidence for density variations in every 
proposed Type cllb supernova for which there are sufficient radio observations to show the 
effect. On the other hand, the Type lb event SN 2008D JSoderberget all l2008h an d the 
Type Ic events SN 2003L JSoderberg etaD l2005h and SN 2009bb (ISoderberg et all l2009h 
show smooth light curves. For SN 2003bg, ISoderberg et al.l (120061 ) estimated that the time 
between the higher densit y outflows was ~ 12( 1000 km s _1 /v w ) yr, taking a timescale of 120 
days to the first feature. iKotak fc Vinkl (120061 ) took a timescale of 150 days and derived a 
period for the mass loss events of ~ 25(200 km s' Vf,,,) yr for SN 2003b g and SN 2001ig. 
The main difference between these estimates is that IKotak & Vinkl (120061 ) used an optically 
determined maximum ejecta velocity o f 15, 000 km s~ l a t 14 da ys to determine the expansion 
of the radio emitting region, whereas ISoderberg et al.l (120061 ) used the SSA model for the 
radio emission to determine a mean velocity of about 40,000 km s^ 1 . As discussed above 
and shown by the VLBI observations of SN 2008ax, the highest velocity ejecta may not be 
detectable at optical wavelengths. The variability timescales for SN 2008ax and SN 2008bo 
are not well defined, but appear to be shorter than in the first two cases. 



As noted by IKotak fc Vinkl (120061 ) , the timescale of decades is consistent with the vari- 
ability timescale of luminous blue variables (LBV) undergoing S Doradus-type variations 
and they propose such stars as progenitors. An example of such a star is AG Car, which ha s 
been followed from its minimum (hotter) state through a cooler phase ( IStahl et al.ll200ll ). 



In its minimum state, the mass loss rate is ~ 3 x 10 5 M yr 1 and v w ~ 300 km s 
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so that A* « 10. This is at the high end, but within the range, of wind densities dis- 
cussed here for the Type c llb events. Howe ver, the radius of AG Car in its compact state 
is 50 R Q = 3.5 x 10 12 cm (jstahl et al. 2001 ). This radius is a factor > 10 time s the radius 
deduc ed here for the progenitor of SN 2008ax and the radii of WR stars. Also, IStahl et al. 



(120011 ) estimate a mass loss rate of ~ 1.5 x 10 



M 



© yr 



and v„ 



150 km s 1 in the 



cool state, so that the wind density is ~ 10 times higher than in the hot state. The density 
contrast is higher if the faster wind sweeps the slower wind into a shell. The shell densities 
are higher than appear present in the Ty pe cllb events and are closer to the shell interaction 
observed in SN 2006jc (IFoley et al.l 120071 ). S Doradus stars with lower luminosities than AG 
Car, which are more likely progenitors of the objects discussed here, have even larger radii. 



Kotak fc Vinkl (120061 ) note that a WR progenitor is possible if the star becomes a such 



a star just before the supernova so it still interacted with mass loss from the LBV phase. 
However this would require the transition < 25 yr before the explosion. In addition, the fast 
WR wind would modify the close in circumstellar medium. Finally, the eje cta masses for SN e 
cllb are not particularly large: estimates ar e 2.5 — 5 M for SN 1996cb (IDeng et al.l 12001 
4 M Q for SN 2003bg jMazzali et alJboOflh . and ~ 3.8 M Q for SN 2008ax jTsvetkov et al. 



20091 ). Unless a considerable amount of matter goes into a black hole central remnant, these 
masses are less than expected for the cores of LBV stars. Overall, the evidence points to 
progenitors for Type cllb supernovae that are less massive and more compact than LBVs, 
and are more consistent with Wolf-Rayet stars. Only a weak variation in wind properties 
is necessary to produce the observed radio variations. However, the possibility that the 
mechanism is similar to that producing S Doradus variations in LBVs cannot by ruled out. 



DISCUSSION 



One expectation of the two types of lib is that the Type cllb events should have a 
lower H mass e nvelo pe so that a large radial extent cannot be sustained. For SN 1993J, 
Wooslev et all (|l994h estimate a H mass of 0.2 ± 0.05 M . For SN 2003 bg, Mazzali et al. 



( 20091 ) estimate 0.05 M of H, in keeping with expectations. However, Qiu et al. ( 19991 ) 
suggest that SN 1996cb has more H than SN 1993J, although IDeng et al.l (1200 ll ) estimate 
0.1 — 0.2 M & of H if the explosion energy is 10 51 ergs. The implied dividing line between the 
types is ~ 0.1 M . 

For lower amounts of H, we expect that the Type cllb objects merge into the Type 
Ib's. A recent example of a Type lb with H is SN 2007Y which showed high velocity 
(~ 10, 000 km s _1 ) H absorption near maximum l ight (IStritzinger et al.l 120091) . The finding 



of H in the spectra of Type lb events is common (jBranch et al.ll2002f ) 



Branch et al.l (12002 ) 
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estimate that a H mass ~ 0.01 M & is typically needed to produce the required optical depth 
in Ha. There may not be a clear separation of the Type cllb and Type lb events, but a 
gradual transition depending on remaining H mass. 

The Cas A supernov a was recently shown to be of Type lib by the spectrum of its light 



echo (IKrause et al.l 120081 ) . which raises the question of whether the proge nitor was compact 



or ext ended. Studies of the supernova remnant have not been conclusive; IChevalier fc Oishi 



(120031 ) argued that the progenitor was extended with the dense wind extendi ng to the stellar 



surfa ce, but there have been arguments for a blue phase before the supernova (IHwang fc Laming 



20091 ). The Type lib identification was made on the basis of the close correspondence of the 
echo spectrum with the time integrated spectrum of SN 1993J. The Ha line in the echo 
spectrum appears t o be closer to the maximum light Ha in SN 1993J than in SN 2003bg 



(IHamuy et al.l 120091 ). but a more detailed comparison with the time integrated spectrum of 



a Type cllb event would be useful. 

The existing observations that allow a distinction between Type cllb and ellb super- 
novae, primarily at radio wavelengths, indicate a clear separation between the 2 types of 
events, but the numbers are small. Future time domain surveys should increase the discov- 
ery rate of Type lib events so that this issue can be resolved. In addition, the early discovery 
of supernovae (< 2 days from explosion) will allow the classification of the events from their 
shock breakout and cooling envelope emission. 

We are grateful to the referee for a helpful report, and acknowledge support from NSF 
grant AST-0807727 (RAC), a Hubble Fellowship (AMS), and Swift GI program NNX09AR05G. 
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Fig. 1. — Peak spectral radio luminosity vs. the product of the time of the peak and 
the frequency of the measurement for supernovae of Types Ib/c and lib. The observed 
supernovae are designated b y the last two dig i ts of the year and let ter (s). References to the 
recent radio obse rvations are poderberg et al.l (120081 ) for SN 2008D. lRoming et al.l (120091 ) for 
SN 2008ax. and ISoderberg et al.l (120091 ) for SN 2009bb; other references can be found in 



Chevalier fc Franssonl ( 120061 ). The dashed lines show the mean velocity of the radio shell if 
synchrotron sel f-absorption is r esponsible for the flux peak; a particle spectral index p = 3 



is assumed (see IChevalierl Il998l ) . 



